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We present a detailed study of the charge transport, optical reflectivity, and thermal transport 
properties of n-type PbSe crystals. A strong scattering, mobility-limiting mechanism was revealed 
to be at play at temperatures above 500 K. The mechanism is indicative of complex electron-phonon 
interactions that cannot be explained by conventional acoustical phonon scattering alone. We ap- 
plied the first order non-parabolicity approximation to extract the density of states effective mass as 
a function of doping both at room temperature and at 700 K. The results are compared to those of a 
parabolic band model and in the light of doping dependent studies of the infrared optical reflectivity. 
The thermal conductivity behavior as a function of temperature shows strong deviation from the 
expected Debye-Peierls high temperature behavior (umklapp dominated) indicating an additional 
heat carrying channel, which we associate with optical phonon excitations. The correlation of the 
thermal conductivity observations to the high temperature carrier mobility behavior is discussed. 
The thermoelectric figure of merit exhibits a promising value of ~ 0.8 at 700K at 



I. INTRODUCTION 

The cubic semiconductor PbSe (space group Fm-3m, 
a=6.125 A) has been attracting scientific attention for 
more that seven decades mainly because of its opti- 
cal properties. For example, the infra-red behavior of 
PbSe has been appreciated since the 1940si and even- 
tually found applications in photodetectors and thermal 
imaging. 2 More recently, the advent of nanoscience raised 
interest in photovoltaic applications based on excitonic 
effects in PbSe nanocrystals£~— Hence, the compelling 
majority of published work on the system involves the 
study of thin films and other nanostructures mainly at 
room temperature and lower. 

Lately, however, theoretical and experimental reports 
have pointed out the appealing characteristics of PbSe 
for higher temperature thermoelectric applications^— 
For example, PbSe melts at a relatively high tempera- 
ture (1080 °C), is composed of earth abundant elements, 
is easy to scale up, and performs better than PbTe at 
900 K, one of the choice thermoelectric materials in the 



temperature regime 600-900 Kj2t— It is noteworthy, that 
all of the above reports (refs 6-10) have highlighted the 
lack of detailed experimental studies on the charge and 
thermal transport as well as basic electronic band struc- 
ture parameters of PbSe, with different dopants and as a 
function of doping level especially at high temperatures. 

Hirahara et al, for example, studied the mobility of 
both n and p type PbSe hot pressed samples up to 
773 K and doping levels well below 10 19 cm~ 3 taking 
into account impurity scattering^ Later, Schlichting and 
Gobrecht repeated the mobility study up to 800 K on 
melt grown crystals of n and p type PbSe extending the 
doping levels up to 3 x 10 19 cm -3 , and concluded that 
electron phonon interactions are dominant/i 3 - Scattering 
was also studied in the framework of defect formation 
through measurements of Hall effect, electrical conduc- 
tivity, and thermopower by Gurieva et al. but the study 
was not conclusive as to which type of defects (e.g. in- 
terstitial, Frenkel etc) dominate.— Alekseeva et al. have 
studied the high temperature properties of p-type PbSe 
and samples with isovalent Cd and Mn ion substitution 
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respectively. 15 Finally, limited attention has been given 
to the valence band structure of PbSe and its effect on 
the high temperature thermoelectric properties . 16 ' 17 

Given the increasing interest in the high temperature 
properties of PbSe and the aforementioned conflicting 
conclusions on the charge transport properties of PbSe, it 
becomes clear that a reliable and consistent experimental 
body of results has to be established. Hence, we initiated 
a study of n-type samples doped with CI exhibiting an 
electron carrier density in the range 7.5 x 10 18 < n < 
3.8 x 10 19 cm -3 . The choice of CI as a dopant is justified 
by its substitutionary action on the Se sublattice leav- 
ing undisturbed the conduction band of PbSe that con- 
sists primarily of Pb p-orbitals*^ 7 - We report electrical 
conductivity, Hall coefficient, thermoelectric power and 
thermal conductivity as a function of doping and tem- 
perature. Furthermore, the optical reflectivity as a func- 
tion of doping at room temperature was studied within 
the framework of a Krammers-Kroning analysis. We em- 
ployed both a parabolic and a non-parabolic (Kane-type) 
band model to extract basic parameters such as the ef- 
fective mass and the Lorenz number at different temper- 
atures and as a function of the doping level. Finally, 
a thorough investigation of the lattice thermal conduc- 
tivity is presented where in addition to acoustical and 
three phonon processes optical phonon contributions are 
required to better account for the observed high tempera- 
ture behavior. We find that the maximum thermoelectric 
figure of merit, ZT, achieved at 700 K is ~0.8 for a carrier 
density of ~ 1.5 x 10 19 cm~ 3 ^ 



II. EXPERIMENTAL DETAILS 

The PbSe crystals were grown by the Bridgman tech- 
nique inside sealed and evacuated quartz ampules with 
one end tapered. The ampules were loaded with high pu- 
rity Pb (99.999%, American Elements), Se (99.999%, 5N 
Plus) and PbCl 2 (99.9999%, Aldrich). Initially the load 
was suspended in the hot zone of the furnace at 600 °C 
for 72 h. Then the hot zone was heated to 1150 °C, the 
load was raised at a higher position outside the hot zone 
and dipped at a speed of ~1.2 mm/h. After growth, 
the ingots were sliced to 8 mm diameter disks with a 
waferizing saw. Subsequently, two of the disks were fur- 
ther processes with a polisher to form a bar of typical 
dimensions ^7x3.5x2.5 mm and to a disk of thickness 
^2 mm. Optical examination of the specimens revealed a 
polycrystalline texture consisting of large single crystals 
oriented at different directions. 

The bar and disk specimens of each growth were sepa- 
rately examined at room temperature for consistency in 
doping by combining Hall effect and a spatial scanning 
Seebeck apparatus. Specimens with large inhomogeneous 
areas of thermopower, pertaining to large carrier density 
variations, were not considered further. Figure 1 presents 
a typical example of a homogeneous and an inhomoge- 
neous pair of specimens. Generally, we have observed 



TABLE I: Specimens of the present study tabulated based on 
PbCl2 content and corresponding carrier, plasmon frequency, 
and mass density. 



ID PbCl 2 %mol n (10 ly cm- a ) uj p (cm" 1 ) p (g/cm 3 ) 



A 


0.05 


0.75 


697 


7.90 


B 


0.10 


0.88 


710 


7.98 


C 


0.20 


1.50 


729 


8.02 


D 


0.30 


3.20 


970 


7.90 


E 


0.40 


3.80 


990 


7.99 


F 




<0.2 




8.09 



that doping levels above 0.4 % mol, corresponding to an 
average Hall concentration, n, higher than 4xl0 19 cm~ 3 , 
produce such high doping inhomogeneities. Therefore, 
the present study was limited to a set of five samples ex- 
hibiting high homogeneity and an average carrier density 
not higher than 3.8xl0 19 cm -3 , see Table I. The mea- 
sured n increases monotonically as a function of increas- 
ing nominal CI concentration. A one to one correspon- 
dence of n versus CI concentration was observed to a rea- 
sonable degree at n < 1.5 x 10 19 cm~ 3 . At high PbCl2 
molar concentrations there is a deviation indicating dop- 
ing action inability of CI in the PbSe lattice. 

The high temperature Hall coefficient was measured 
in a home-made high temperature apparatus, which pro- 
vides a working range from 300 K to 873 K. The sam- 
ples were press-mounted and protected with argon gas to 
avoid possible oxidization at high temperature. The Hall 
resistance was monitored with a Linear Research AC Re- 
sistance Bridge (LR-700), and the data were taken in a 
field of ± 1 T provided by an Oxford Superconducting 
air-bore magnet. 

The electrical conductivity, a, and Seebeck coefficient, 
S, were measured simultaneously on the bar shaped spec- 
imens in a ULVAC-RIKO ZEM-3 system. The specimens 
were protected in a helium atmosphere (~ 0.1 atm) while 
the furnace of the instrument was cycled from room tem- 
perature to ~ 700 K and back. No thermal hysteresis was 
observed with thermal cycling. 

The disk-shaped specimens were used to determine the 
thermal diffusivity as a function of temperature and dop- 
ing in a NETZSCH LFA 457 Microflash instrument. Sub- 
sequently the thermal conductivity, k, was estimated by 
the relation k = DC p p, where D is the thermal diffusiv- 
ity, C p is the heat capacity under constant pressure and 
p is the mass density of the specimens. C p was approxi- 
mated by the formula 0.171 + (2.65 x 10~ 5 )T. 9 All charge 
and thermal transport measurements were performed in 
the same specimen direction. 

Room temperature infrared reflectivity (IR) measure- 
ments were performed on finely-polished PbSe samples 
using a Bruker 113V FTIR spectrometer. The spectra 
were collected in the 100-3000 cm -1 spectral region with 
a resolution of 2 cm -1 at nearly normal incidence. The 
reflection coefficient was determined by a typical samplc- 
in-sample-out method with a mirror as the reference. The 
Im(e) and Im{— 1/e) spectra (where e is the complex di- 
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FIG. 1: Scanning Seebeck results on two coins doped with 0.3% mol PbCb (a) and 0.5% mol PbCbj (b). The small inhomo- 
geneous regions randomly observed in (a), which are acceptable for an ingot material, are extended and create steep gradients 
with increasing PbCbj concentration (b). 



electric function) were derived from the Kramers-Kronig 
transformation^ 



III. 



RESULTS AND DISCUSSION 



A. Charge Transport Measurements 

The electrical conductivity, a of samples A-E is de- 
picted in Fig. 2. The values of a are increasing with 
increasing doping, i.e. moving from specimen A to E, 
at any given temperature. At room temperature a val- 
ues as high as 3500 S/cm can be realized for a doping 
level of 3.8 xlO 19 cm -3 , indicative of relatively high mo- 
bilities, \i. For all specimens a monotonic decrease in a 
with increasing temperature is observed. Since a ~ n/ie, 
the functional dependence of a in temperature may re- 
sult from the temperature dependence either of n or from 
factors limiting /i. 

To elucidate the behavior of a we performed temper- 
ature dependent Hall effect studies. The Hall coeffi- 
cient, Rh, for samples A, C, and E, is plotted in Fig. 
3a as a function of temperature. Rh is almost tem- 
perature independent. Considering Rh — 1/ne it fol- 
lows that the specimens retain the same carrier den- 
sity up to 700 K. A more accurate description of Rh 
should take into account the non-parabolicity of the 
bands and the statistical anisotropy due to scattering. 
In such a case Rh = Ar/ne^ where r is the statisti- 
cal anisotropy factor that varies from 1 to 37r/8, and A 
is the energy surface anisotropy factor that is equal to 
A = 3K{K + 2)/ \2K + 1) 2 M For PbSe, K = 1.75 2 =± and 
therefore the product Ar varies from 0.97 to 1.14 depend- 
ing on the value of r. For simplicity we have kept Ar = 1, 
since fundamentally the conclusions of our study are not 
distorted by such an assumption. 

Figure 3b, presents the Hall mobility, fin ~ Rh&, as 
a function of temperature for specimens A, C, and E. 
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FIG. 2: Electrical conductivity as a function of temperature 
for samples A-E (see Table I for carrier concentrations). The 
conductivity increases with increasing doping at any given 
temperature and no irreversible effects are observed with ther- 
mal cycling. 



(Specimens B and D exhibit similar behavior but are not 
included in the discussion for clarity). The room tem- 
perature [in values are quite high and drop from ~1200 
cm 2 /Vs for specimen A to ~600 cm 2 /Vs for specimen 
E. The decrease with doping reflects increasing e-e in- 
teractions with band filling. The high values of fj,n are 
not retained at high temperatures. A rapid decrease is 
observed that limits fi H to —100 cm 2 /Vs at 700 K. The 
rapid decrease of /j,h is clearer on a logfin — l°gT plot 
where two linear regions show up, see Fig. 4. In the 
temperature regime 300 < T < 480 K the slopes of the 
curves are ~-2. However, for T >520K the slopes in- 
crease to —-3. Table II, summarizes the fitting results. 
With increasing doping concentration the slopes assume 
lower values. 
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FIG. 3: a) Hall coefficient as a function of temperature for 
specimens A, C and E. Note that the Hall coefficient is almost 
temperature independent, b) Hall mobilities for specimens 
A, C, and E as a function of temperature. Despite the high 
room temperature values a drastic decrease is observed at 
high temperatures. 



TABLE II: Characteristic slopes of the log/in -logT for sam- 
ples A, C, and E for two different temperature regimes. 

300> T >480 K 550> T >750 K 
ACE ACE 
-2.2 -2.1 -1.75 -3.2 -3.2 -2.8 



Electronic scattering in semiconductors gives rise to 
distinct temperature dependencies of the carrier mobil- 
ity. Typical processes considered include: (i) electron- 
phonon scattering due to thermal vibrations of the lat- 
tice, (ii) ionized impurity scattering, (hi) scattering from 
high frequency vibrations of the lattice (optical phonons), 
and (iv) scattering by neutral impurities. On the as- 
sumption of parabolic bands, thermal lattice scattering 
causes fin to scale as (j,^ ~ T 3 / 2 M- Non parabolicity 
of the bands, as in the case of Si,— generally increases 
the value of the temperature exponent to ^2-2. 5. 24 In a 
similar fashion the low temperature exponent in Fig. 4 



FIG. 4: a) log/in ~ logT plot showing two regions of linear- 
ity one below 480 K and the other above 520 K. The high 
temperature region is characterized by a much steeper slope 
reflecting a strong electronic scattering mechanism at play, b) 
Mobility as a function of temperature scaled according to eq. 
1, see text. A linear region is observed that extends from 500 
to 700 K. 



can be attributed to acoustical phonon scattering of elec- 
trons, in agreement with previous results on highly doped 
PbSe samples.— The second mechanism, i.e. ionized im- 
purity scattering, constitutes a positive contribution to 
mobility 2 ^ and is approximated as [in ~ 7 13 / 2 £k Ion- 
ized impurity scattering shows up mostly at low tempera- 
tures since the increasing thermal velocity of carriers with 
increasing temperature effectively screens the Coulomb 
potential of impurities^ Finally, mechanism (iv) con- 
tributes only weakly to scattering and does not have a 
significant temperature dependence^ Therefore, in or- 
der to elucidate the behavior of the mobility at T > 500 
K we have to consider more complex processes such as 
polar optical phonon scattering. 

The temperature dependence of the electron mobility 
due to scattering by high frequency optical phonons of 
the lattice follows the analysis of Fortini et al.^ that 
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showed a temperature dependence of the form: 

/i « CT^V - 1)GW (1) 

where, C is a constant of proportionality that involves 
the static and high frequency dielectric constants and 
the frequency of the longitudinal optical phonons, r\ is 
the reduced chemical potential and G[rf) is a function 
that assumes values from 0.65 to 1.7 depending on the 
excitation energy of the optical phonons. We point out 
that Mott and Gurney reached a surprisingly similar re- 
sult, i.e. a formula containing an activation energy, from 
a completely different standpoint in the case of ionic po- 
lar crystals: 27 



fj. ~ (e 



e/T 



1) 



(2) 



where, is a characteristic temperature between 300 to 
800 K. In Fig. 4b we present high temperature mobility 
data as a function of T 1 / 2 (e T °l T — 1), where we define To 
as a characteristic activation temperature corresponding 
to an energy of ~50 meV. The linear region observed 
in this case indicates that complicated scattering mecha- 
nisms, possibly involving optical phonons, are at play at 
high temperatures. 

However, the data are not conclusive as to weather an 
electron-optical phonon scattering is taking place. This is 
because in PbSe, where the lowest conduction band min- 
imum is at the L point of the first Brillouin zone ([111/, 
direction)," and hence is highly degenerate, phonons may 
scatter electrons transferring them from one valley to 
another. 25 This intervalley scattering mechanism can in- 
volve both acoustical and optical phonons and leads to 
a temperature dependence with similar characteristics as 
in eq. 1 and 2. 25 However, in this case many more pa- 
rameters are involved such as the intervalley energy sep- 
aration, the intervalley deformation potential etc.j=s and 
hence a theoretical study is required to reach a definite 
conclusion as to the exact nature of the mobility-limiting 
mechanism at high temperatures. In any case, the data 
clearly point to deviations from simple electron phonon 
interactions at temperatures above 500 K and doping lev- 
els n > 7 x 10 18 cm -3 . Here we note that Schlichting et 
al. presented measurements of highly doped n > 3 x 10 18 
cm -3 PbSe specimens only up to 500 K 13 which may 
have led to the erroneous conclusion for a simple domi- 
nant acoustical phonon scattering mechanism regardless 
of the doping level. At lower n and T > 500K mobility 
data are plagued by bipolar diffusion which masks other 
interactions^ 

The Seebeck coefficient, S for specimens A-E is pre- 
sented in Fig. 5 as a function of temperature. For all 
samples, i.e. for any doping level, S is almost linearly 
decreasing from 300 to 700 K assuming higher absolute 
values with increasing temperature. At 300 K the abso- 
lute value of S is decreasing from sample A to E consis- 
tent with increasing carrier density. The same behavior 
is observed at all temperatures. 

Assuming parabolic bands and an energy independent 
relaxation time then at any temperature S(n) is de- 
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FIG. 5: Seebeck coefficient as a function of temperature for 
samples A-E. 



scribed by a unique effective mass value, m* , and for any 
n the S(T) depends on the temperature dependence of 
m* i 28 ' 29 The aforementioned assumptions lead to simple 
formulas that have been explained in detail elsewhere^ 
and can be used to construct Pisarenko plots, i.e. S(n) 
diagrams at different temperatures. Such plots are de- 
picted in Fig. 6 for three different temperatures, 300 K, 
500 K, and 700 K and the corresponding m* values 0.28, 
0.35, and 0.41 that allowed fitting of the data (solid lines). 
The effective mass value is increasing as a result of an 
increasing band gap, E g , with increasing temperature^ 
consistent with the results of k ■ J? theory.— The tem- 
perature dependence of m* is plotted in Fig. 7a and 
compared to E g {T). For PbSe E g (300K) « 0.275 eV 9 
and dEg/dT « 4 x 10~ 4 eV/K. 31 It is evident that in the 
crude approximation of the parabolic band model m* fol- 
lows very close the rate of the band gap increase. 

As mentioned earlier, the above values are only crude 
approximations since the electronic band structure of 
PbSe close to the Fermi level is non parabolic. 7 Non 
parabolicity pertains to a non-spherical Fermi surface 
shape, and hence the dispersion relation depends on the 
spatial direction. In general, a non parabolic energy band 
dispersion can be expanded in a power series^ 



q=2 



(3) 



where the coefficients \ q are defined from the following 
relation^ 



_ h 2 1 d q k 2 
q ~ 2^^! ( ^" ) 



(4) 



The first order non-parabolicity approximation ignores 
all terms with q > 3 and thus eq. 1 is greatly simplified 
to: 



, 9 2m , , , 
k 2 = — e(l + Ae) 



(5) 
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FIG. 6: Pisarenko plots at 300 K, 500 K, and 700 K. The solid 
lines represent fittings of the data points on the assumption 
of a parabolic conduction band and an energy independent 
scattering time. The only adjustable parameter in the cal- 
culation was m*. Note that the higher the temperature the 
higher the effective mass. 



In this approximation the coefficient A is usually taken 
to be equal to the inverse of the band gap, E g ££ Conse- 
quently, all galvanomagnetic coefficients can be expressed 
as functions of the generalized Fermi integrals % L\ defined 



by the equation>si 

i L\{r}^)= r(~^)z\z + (3z 2 y(l + 2/3z) l dz (6) 
Jo az 

where fo(rj,T) is the Fermi distribution function, 77 is 
the reduced chemical potential, z = e/k B T, and /3 = 
Xk B T = k B T/E g . 

Here, we are especially concerned with the expression 
for S. Fitting S yields rj values that can be utilized in 
the calculation of other parameters in the same fashion 
as with assuming a parabolic band.— S is defined by i 21 ' 32 



l LL 2 -r,(?LU) 



"L 1 



(7) 



Using the r\ values extracted from fitting the S(n) data 
with eq. 7 and in conjunction with the expression for the 
carrier density (with a unity Hall factor)^ 



1 

3^ 



2m*k B T 



WO 



(8) 



the density of states effective mass, m?, can be calculated 
at any given pair of n and T values. The dependence 
of m* d on n is another substantial difference from the 
parabolic band model and is supported experimentally 
in lead chalcogenides i^2r— 

The results of the fitting process using eq. 7 and 8 are 
depicted graphically in Fig. 7b for 300 K and 700 K. The 
extracted values are lower compared to those calculated 
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FIG. 7: a) Temperature dependence of the effective mass as 
calculated from the parabolic band model and comparison 
with the temperature dependence of the band gap. Note that 
the results are independent of the carrier density, b) The 
effective mass calculated from Seebeck coefficient data in the 
first order non-parabolic approximation. 



from the parabolic band model, see Fig. 7a, but closer to 
0.21, the textbook value for PbSe^ Both at 300 and 700 
K m* d is monotonically increasing with increasing doping 
as discussed above. 



B. Optical Reflectivity: independent measurement 
of the effective mass 

The study of the reflectivity as a function of energy in 
the infrared part of the electromagnetic spectrum yields 
useful information about basic materials parameters that 
are related to the electronic band structure, such as the 
effective massJ^ This is because the reflectivity, R, de- 
pends on the contribution of both the bound and the 
free electrons to the real part of the complex index of 
refraction. Therefore, we have performed infrared reflec- 
tivity measurements in our samples as an independent, 
yet direct, method of probing the effective mass and thus 
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FIG. 9: The density of states effective mass as a function 
of increasing carrier density extracted from: thermopower 
measurements analyzed with a single parabolic band model 
(magenta line), thermopower measurements analyzed in the 
first order non-parabolicity approximation (* marks), op- 
tical reflectivity measurements with the assumption that 
e 3000cm -i ~ too (solid green squares), optical reflectivity mea- 
surements with a correction of ~ 25 for all specimens (blue 
x marks). 



FIG. 8: Room temperature reflectivity spectra as a function 
increasing frequency. Spectra were shifted vertically by 0.4 for 
clarity. The reflectivity minimum shifts to higher values from 
sample A to E as a result of increasing carrier concentration. 



compare the results with those of the transport models 
as discussed in the previous section. 

The room temperature infra red optical reflectivity as 
a function of incident radiation wavelength is depicted in 
Fig. 8 for samples A through E. It is readily seen that 
the minimum, associated with the plasma frequency, up 
and hence the carrier concentration, is steadily increasing 
from specimen A to E reflecting the increasing doping. 
Accordingly, the Kramers-Kronig transformation yielded 
blue-shifted absorption peaks in 7m(l/e), in the region 
>700 cmX, that were used to accurately determine the 
plasma frequency, up. The transverse optical - longitu- 
dinal optical splitting in the reflectivity spectra of PbSe 
occurs in the 34-114 cm -1 region^ The latter makes 
plasmon-phonon effects negligible in the present study. 
Therefore, the reflectivity minima in Fig. 8 are mainly 
determined by the contribution of free carriers. 

The plasma frequency is related to basic materials pa- 
rameters, such as the electric susceptibility (or optical) 
effective mass, ? 



through the relation: 



(9) 



where n is the carrier concentration, f M is the high fre- 
quency dielectric constant (a measure of the bound elec- 
tron contributions to the dielectric function), eo the vac- 



uum permeability, and e is the electron charge. The 
expression for the dependence of R on frequency, u), at 
nearly normal incidence: 



i?H = ( 



(10) 



where e(u>) is the complex dielectric function. Using the 
experimental R values at 3000 cm -1 and eq. 10 we cal- 
culated the value of (eoo « tzooocm- 1 )- We have ob- 
served that too is decreasing from 25 to 17 with increasing 
n. However, there is no physical reason for a changing 
contribution of the bound electrons with n in PbSe, and 
thus, we have taken eoo ~ 25 for all samples. 

The values of up (see Table I) were derived from the 
peak value of Im(—l/e) obtained by Kramers-Kronig 
transformation. Using eq. 9, with the Hall-cffcct ex- 
tracted n we have evaluated m* p . The calculated val- 
ues of m* p can be transformed to m* d by using the 
relation ; 21 ! 35 



1 + 2K 



(11) 



where N m is the number of equivallent conduction band 
ellipsoids in the first Brillouin zone and K is the ellip- 
soid anisotropy factor which for PbSe takes the value 
K = 1.75i2i Our results are plotted in Fig. 9 as the 
green squares and compared to the extracted m* d values 
from thermopower data analysis (* marks). We observe 
that the optical measurements support a strongly in- 
creasing m* d with increasing n in agreement with the non- 
parabolic nature of the conduction band of PbSe. The 
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FIG. 10: Total thermal conductivity as a function of temper- 
ature. 



agreement between the optically extracted and charge 
transport extracted m* d values is excellent for eoo ~ 25 
(blue x marks) . Obviously, the divergence of the data at 
high n are associated with an underestimation of the 
values with heavy doping. 



C. Thermal Transport 

The total thermal conductivity (k = n e + Ki a t, where 
Ki a t is the lattice part and K e the free carrier contribu- 
tion) as a function of temperature is presented in Fig. 10. 
The room temperature values start from ~2.4 W/mK 
for sample A and increase with increasing doping to ^4 
W/mK for sample E. Rising temperature increases the 
electron-phonon and phonon-phonon interactions which 
cause k to decrease. Interestingly, k remains above 1 
W/mK at all temperatures, despite the strong mobility- 
reducing mechanism that is in effect above 500 K and 
should also limit the heat carrying efficiency of carriers. 
Therefore, we conclude that the latter may result from 
an increased lattice contribution at high temperatures, 
i.e. another mechanism contributing to heat conduction. 

Kiat is estimated indirectly by use of the Wiedemann- 
Franz relation, K e — LaT, where L is the Lorenz number 
and subtraction from the total. The temperature depen- 
dence of L is critical in the proper calculation of n e . A 
good approximation that yields reasonable results is the 
assumption of a parabolic band, electron acoustic phonon 
interaction only and a constant relaxation time^^ Fig- 
ure 9a presents the results of such a calculation for sam- 
ples A, C, and E. The calculated Lorenz number values 
are below the metallic limit (L « 2.45 x 10~ 8 WQ/K 2 ) 
and decrease with decreasing doping and increasing tem- 
perature. 

In the first-order non parabolic approximation the 



) 



(12) 



where l Lj are the integrals defined by eq. 6 and are 
functions of rj. The temperature dependence of rj is ex- 
tracted by fitting S(T) data, see Fig. 5, with eq. 7. 
The results are plotted as a function of temperature in 
Fig. lib. Comparing with the results of the parabolic 
model, it is evident that the first order non-parabolicity 
leads to higher Lorenz number values for the same doping 
level. At high doping (sample E) the room temperature 
value is slightly larger than that of the metallic limit, 
L « 2.45 x l0~ a Wfl/K 2 . This may be due to an inade- 
quacy of the non-parabolicity approximation (A ~ E^ 1 ) 
at high doping. Both models, however, exhibit essentially 
the same functional dependence with respect to temper- 
ature and therefore the temperature dependence of Ki a t 
remains the same. 

Due to the overestimation of the Lorenz number of 
the non-parabolic model at high doping we have used 
the predictions of the parabolic model to extract lattice 
thermal conductivity as a function of temperature. The 
results are depicted graphically in Fig. 12a. For clarity 
we use the results for samples A and E only. Curves for 
samples B-D have similar values and temperature depen- 
dence, leading us to conclude that Ki a t is independent of 
doping consistent with the low concentration of PbCl2 
used in this study. 

Generally the temperature behavior of Ki a t when umk- 
lapp processes are dominant is K; at ~ T^ 1 . This is the 
case of PbTe4£ In the case of PbSe, however, we find 
that Ki a t T 1 - 5 with 5 « 0.2 - 0.23. In Fig. 12b ni at is 
presented as a function of lOOO/T 8 . The solid line is a 
linear fit of the data confirming the scaling behavior. In 
order to exclude such a behavior stemming from fitting 
artifacts (Lorenz number calculation) or from doping, de- 
spite the dilute CI concentration, we grew separately a 
pure, undoped PbSe single crystal, here referred to as 
sample F in Table I. (The same crystal was also used 
in a previous study^). The carrier concentration in the 
undoped crystal was determined through Hall measure- 
ments to be < 2 x 10 18 cto~ 3 . In such a case the contribu- 
tion of free carriers is minimized. In confirmation of the 
scaling behavior mentioned above the Ki a t of specimen F 
exhibits a similar temperature dependence with 8 ks 0.23 
and a room temperature value of ~1.9 W/mK. 

Usually a positive S is attributed to optical phonon ex- 
citations that provide an extra heat carrying path42r— 
This raises the thermal conductivity at high tempera- 
tures compared to the simple 6 = behavior. Steigmcicr 
and Kudman, used the well known result of three- phonon 
processes at high temperatures: 



Kiat 



1 h ' -f 2 T 



(13) 



in combination with temperature dependent «/ at accu- 
rate data of several III-V compounds received with the 
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FIG. 11: a) Lorenz number for samples A, C, and E calculated 
as a function of temperature on the assumption of a parabolic 
conduction band, electron-acoustical phonon interactions only 
and constant relaxation time, b) Same as in a) however the 
first order non-parabolicity approximation is used, see text. 



FIG. 12: a) Lattice thermal conductivity extracted using the 
L P ar, see Fig. 11, temperature dependence. Notice that Kiat 
is doping independent, b) Scaling of the lattice thermal con- 
ductivity as ~ T 0,8 . The dashed line is a linear fit (sample 
A), used here as a guide to the eye. 



flash diffusivity - heat capacity method.- 3 They con- 
cluded that eq. 13 (where, ks,h,a, p,0D,^f are the Boltz- 
mann constant, Planck's constant, the lattice parameter, 
the density, the Debye temperature, and the Griineisen 
parameter respectively) is valid only when 7 is tem- 
perature dependent, indicating the presence of optical 
phonons contributing to heat conduction. A similar anal- 
ysis on our PbSe specimens leads to the same conclusion 
on 7. Interestingly, a more simplified model developed 
by Dugdale et al41 with a similar physical basis as the 
three phonon model, was previously found adequate to 
describe the temperature dependence of Ki a t on the as- 
sumption of a temperature dependent 7^ 

To separate the contribution of acoustical phonons 
from the Ki a t we have calculated the phonon thermal 
conductivity by assuming a phonon scattering relaxation 
time t(x)^L 

t(x) 1 = t d 1 + Tp 1 = Au A + CTw 2 (14) 

where we have considered only contributions from point 



defects (D index) and umklapp processes (P index). In 
the above u is the phonon frequency. Hence the acousti- 
cal phonon contribution to K/ Qt takes the form: 

k B faT 3 f e °/ T x 2 e x 
Kl ' ac ~ 2tt 2 v s ( h ' J Dx* + E> (e* - l) 2 

where v s the speed of sound, x — Tiw/kBT the dimen- 
sionless variable of the phononic energy, D is a tem- 
perature independent constant and E' — CTik^T jK) 2 . 
Since, optical phonons presumably appear at tempera- 
tures above Qb we have used literature 2 -^ Ki a t values 
at T < 8r> = 170K— , i.e. at temperatures where only 
acoustical phonons dominate, to extract the constant C . 
Subsequently E', and n^ ac where calculated. The results 
are depicted in Fig. 13a along with «; ai of sample F. 9 Ev- 
idently, there is a considerable deviation which increases 
with rising temperature. 

Interestingly, Alekseeva et al4^ have previously ob- 
served similar thermal conductivity scalings on the 
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FIG. 13: a) Temperature dependence of the lattice thermal 
conductivity of PbSe (specimen F, Table I, red triangles) com- 
pared to the expected behavior on the assumption of umklapp 
processes only (dashed blue line), b) The contribution of op- 
tical phonons extracted from the difference of the curves in a) 
(green crosses) and also by fitting the measured PbSe lattice 
data of specimen F to eq. 16 (blue x symbols), see text for 
details. Remarkably, the dependence of ni, op on temperature 
is the same regardless of the model used. 



lighter lead chalcogenides, i.e. PbSe and PbS. The au- 
thors formulated the relationship^ 



K lat — K i.ac T Kl , 



T 



+ l3(e k B T - e k B T ) 



(16) 



where Ki iQp is the contribution of optical phonons to 
the lattice thermal conductivity, a' and /3 are fitting 
constants, and ujq and oj\ are the minimum and maxi- 
mum optical phonon frequencies. Equation 16 was ex- 
tracted based on experimental phonon spectra of lead 
chalcogenides 4£ 

The temperature dependence of the contribution of op- 
tical phonons can be calculated either by the subtraction 
K iat — K i,ac, see eq. 15, or directly fitting Ki at with Alek- 
seeva's model (eq. 16). Using data received on sample 
F we have employed both models. Alekseeva's model re- 



sulted in hu>o ~ 12.4 meV, Tiwi — 89.4 meV, and /3 ~ 0.47 
W/mK. (Note that T in Fig. 4b amounts to ~ 50 meV, 
i.e. the average of Hujq, Hlui minimum and maximum op- 
tical phonon energies as suggested by eq. 16). The con- 
tribution of optical phonons to the total lattice thermal 
conductivity as a function of temperature is graphically 
depicted in Fig. 13b. A comparison of the results gives 
a rough 7 % disagreement between the values predicted 
by the two models at all temperatures, with Alekseeva's 
model lying lower. Nevertheless, a striking similarity in 
the functional form of ni_ op with respect to temperature 
is observed, where a broad peak value appears at ~450 K. 
It is interesting that this temperature correlates with the 
transition temperature of the electronic Hall mobilities 
(see Fig. 4a) progressing from a predominant acousti- 
cal phonon scattering region (300-450 K) to a region of 
stronger, possibly optical phonon, scattering (T>520 K). 
We would like to point out that the peak is not related to 
a maximum contribution of optical phonons. Since the 
contribution from acoustical phonons is rapidly decreas- 
ing (~ T~ 4 ) the percentage contribution of the optical 
phonons to the total thermal conductivity is constantly 
increasing with increasing temperature giving rise to the 
T 1 ^ 5 dependence. 

The physical origin of the optical phonon contributions 
identified here is currently unknown but it may be associ- 
ated with the increasing displacement of Pb atoms from 
the octahedron center in the rock salt structure discov- 
ered recently in PbQ (Q=S, Se, Te)<££ Interestingly, even 
in the heaviest lead chalcogenide, i.e. PbTe, anharmonic 
contributions are present in both the charge and the ther- 
mal transport properties. Feit et al^i identified signifi- 
cant polar optical phonon contributions in n-type PbTe 
and more recently Delaire et ali^ identified a strong cou- 
pling of the transverse optic mode with the acoustical 
longitudinal mode which is believed to overall keep the 
thermal conductivity \ow££ However, in contrast to PbS 
and PbSe the lattice thermal conductivity of PbTe ex- 
hibits a 1/T behavior at high temperatures, i.e. 5 = 0. 45 
The above emphasize that despite their striking struc- 
tural similarity and simplicity, lead chalcogenides present 
challenging electronic structure complexity that varies 
from Te to Se and likely S. In any case despite the in- 
creasing experimental evidence for significant participa- 
tion of interactions involving high frequency out of phase, 
i.e. anharmonic, phonons in the heat conduction process 
in a wide and diverse variety of compound a 10 ' 44 ~ 46 ' 52 i 53 
leading to either an increase or a decrease of Ki a t there is 
still no adequate theoretical treatment. 

The ZTfi£ is depicted as a function of temperature 
in Fig. 14. The maximum value (0.8 at 700 K) is as- 
sumed for sample C (n~ 1.5 x 10 19 cm -3 ). Given the 
high temperature mobility reducing mechanism, the ZT 
values reached here emphasize that PbSe is a promising 
thermoelectric material. We note that chemical substi- 
tutions on the Se sublattice with inexpensive and highly 
abundant S have led to ZTs as high as 1.3 at 900 
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FIG. 14: The thermoelectric figure of merit as a function of 
temperature for all samples. The maximum ZT at 700 K is 
observed ~ 1.5 x 10 19 cm -3 , inset. 



IV. CONCLUDING REMARKS 



The first order non-parabolic model, for carrier concen- 
trations < 1 x 10 19 cm~ 3 and approximating the non- 
linear coefficient with the inverse of a linearly increasing 
band gap, was shown to describe satisfactorily basic elec- 
tronic structure parameters of PbSe such as the density 
of states effective mass. The latter was independently ex- 
tracted by optical reflectivity measurements and found to 
be in good agreement with the charge transport results. 
At high temperatures and higher carrier densities a bet- 
ter non-parabolic approximation of the conduction band 
is necessary. Finally, our results indicate a great poten- 
tial of PbSe for applications as a thermoelectric mate- 
rial at high temperatures. The involvement of optical 
phonons in conducting heat at high temperatures imply 
that in PbSe the lattice thermal conductivity at 700-900 
K is higher in relative terms than in PbTe where opti- 
cal phonons are less important. Therefore, strategies for 
reducing thermal conductivity to raise ZT should take 
these findings into account and be aimed at creating op- 
tical phonon scattering mechanisms. 
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